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A dynamical phase transition from reversible to irreversible behavior occurs when particle suspen- 
sions are subjected to uniform oscillatory shear, even in the Stokes flow limit. We consider a more 
general situation with non-uniform strain (e.g. oscillatory channel flow), which is observed to exhibit 
markedly different dynamics. Self-organization and shear-induced migration only partially explain 
the delayed, simultaneous onset of irreversibility across the channel. The onset of irreversibility is 
accompanied by long-range correlated particle motion. This motion leads to particle activity even 
at the channel center, where the strain is negligible, and prevents the system from evolving into a 
reversible state. 
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I. INTRODUCTION 

Suspensions of non-Brownian particles were recently 
shown to flow irreversibly when subjected to sufficiently 
high, uniform oscillatory strain [TJ even at low Reynolds 
number (Stokes flows). The significance of this phe- 
nomenon is enhanced by the fact that time reversibil- 
ity is one of the most basic symmetries of physics, and 
Stokes flows have long been used to demonstrate hydro- 
dynamic reversibility. Furthermore, a distinct threshold 
strain amplitude was found above which particles fail to 
return to their original positions upon flow reversal. This 
system provides one of the few confirmed examples of a 
non-equilibrium phase transition [2]. Since suspensions 
are an important class of fluids in chemical, biological, 
and industrial applications, the appearance of complex 
dynamics in such flows is significant. These experiments 
have led to a number of theoretical works aimed at ex- 
plaining the non-equilibrium phase transition [3HS]. 

Random, irreversible particle interactions provide a 
mechanism for particle self-organization, which leads to 
a divergent relaxation time at the onset of irreversibil- 
ity. Similarities have been demonstrated between this 
transition and plastic depinning phenomena [6]. Uni- 
formly sheared particle suspensions can also exhibit self- 
organized criticality (SOC) [7] when random interactions 
and gravitational settling compete. 

It has not yet been determined how irreversible parti- 
cle motion develops in suspension flows that are strained 
inhomogeneously, a far more general situation. In this 
paper, we consider the onset of irreversible behavior for 
dense suspensions (40% by volume) in an oscillatory 
channel flow. Here, the strain is largest near the walls 
and negligible at the center, so some regions may locally 
exceed the critical strain amplitude, while others do not. 
We find that the transition to irreversibility is markedly 
different compared to the Couette flow case studied pre- 
viously. Self-organization and shear-induced migration 
contribute to a delayed, simultaneous onset of critical 



behavior across the channel. However, non-local particle 
activity generated in high strain regions near the chan- 
nel walls prevents the system from becoming reversible 
even where the strain is negligible, i.e. at the channel 
centerline. In contrast to the Couette flow case, the fluc- 
tuation relaxation time does not diverge at the transition. 
These experiments provide important guidance for mod- 
eling and suggest additional experimental investigations. 



II. EXPERIMENTAL METHODS 

Polymethylmethacrylate (PMMA) particles with a di- 
ameter d = 220 ± 10 fim are suspended at a bulk average 
volume fraction <f> — 0.4 in an index of refraction and 
density matched solvent with a viscosity /i « 3000 cP 
[H [HI IS]- The large particle size and solvent viscosity 
render the particles essentially non-Brownian. A small 
amount of fluorescent dye (Rhodamine 6G) is added to 
the fluid for visualization purposes. The channel consists 
of a slot (width 2a = 0.5 cm, depth 26 = 1.5 cm and 
length L — 15.0 cm) machined into an optically clear, 
cast acrylic block. The walls are polished to remove tool 
marks and an acrylic plate is bonded to the block to close 
the channel. The suspension is driven in an oscillatory 
manner using a syringe pump with a low viscosity silicon 
oil. The interface between the oil and suspension occurs 
well outside the channel test section [Fig. flja)]. 

A laser sheet (514 nm argon ion laser) is used to probe 
a thin slice (« 100 /im) in the bulk of the suspension. The 
interrogation region (w 1 cm long) is centered along the 
length of the channel far from the ends at mid-depth. 
A CCD camera with a longpass filter (> 550 nm) im- 
ages the fluorescent dye in the interstices of the sus- 
pension [51 [TU] allowing the particles to be tracked [Fig. 
Wb)]. The pump is driven in an approximate square 
wave pattern (see Videol [11]) with a peak flow rate of 
1 ml/min corresponding to a maximum Reynolds num- 
ber of Re < 10 -3 , which ensures Stokes flow conditions. 
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FIG. 2: (color online) (a) Local strain amplitude estimated 
from the measured particle displacement profiles (inset), (b) 
The steady-state particle concentration profiles show signifi- 
cant shear-induced migration at large wall strain amplitudes. 




FIG. 1: (color online) (a) The apparatus uses a syringe pump 
to drive an oscillating suspension flow. A laser sheet probes a 
thin plane, which is imaged from below by a CCD camera with 
a longpass filter, (b) The laser sheet excites fluorescent dye 
in the interstices between particles, making them visible, (c) 
Particle tracks overlaid on a multiple exposure image sampled 
once per period for 35 cycles demonstrates the irreversible 
behavior especially in high strain regions near the wall. 



Prior to each experiment, the suspension is mixed thor- 
oughly with a magnetic stirrer bar to provide consistent 
initial conditions. 

For oscillatory motion, the strain amplitude, which 
governs the mobility of the particles, is varied by adjust- 
ing the volumetric displacement amplitude of the recipro- 
cating pump. The suspension is imaged once per period 
for 500 cycles to capture net deviations from reversible 
particle trajectories measured using particle tracking al- 
gorithms [12] • Afterwards, higher speed imaging (5 fps) 
provides the steady-state concentration profile, 4>(y), and 
the measured displacement profile, S(y), is used to obtain 



the strain amplitude, 7 = dS/dy [Fig. pjj. (All mea- 
surements are averaged about the channel centerline to 
improve the statistical sampling.) The overall strain is 
parameterized by the largest strain in the system, 70, 
which occurs at the wall. 



III. RESULTS 

Irreversible phenomena such as shear-induced migra- 
tion have long been observed in suspensions 13 15 1. Nu- 
merous experiments and simulations have been aimed at 
understanding particle migration in pressure-driven pipe 
and channel flows [T6l - fT8] . Two regimes of migration are 
observed in oscillatory pressure-driven systems |19H21j . 
where large oscillation amplitudes cause migration to the 
centerline, and small amplitudes produce a slight migra- 
tion toward the wall. In our study, significant migration 
toward the centerline is observed at large oscillation am- 
plitude, with negligible migration for the smallest strain 
amplitudes tested [Fig. [2^b)]. While the bulk effects 
of these phenomena are well documented, their conse- 
quences for particle self-organization and the onset of ir- 
reversibility are unknown. 
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FIG. 3: (color online) (a) For fixed wall strain amplitude, the 
particle diffusivity D x is largest near the channel walls, and 
also rises with the wall strain amplitude for any fixed position, 
(b) A subset of the same data shows that the development of 
irreversibility is simultaneous throughout the channel. 



A. Steady-state Fluctuations 

The mean square particle displacements (MSDs) mea- 
sure the degree of irreversibility or particle "activity". 
We use the stream-wise diffusivity, D x — (Ax 2 )/2N, 
(MSD per cycle), and we omit transients occurring in 
the first ~ 200 cycles when considering steady-state mea- 
surements. For small wall strain amplitude, the particle 
trajectory deviations are negligible and thus reversible. 
As the wall strain increases, we observe a rapid rise of 
the particle diffusivity in the high strain regions near 
the walls, and a markedly slower increase near the cen- 
ter (Dy is somewhat smaller, but shows similar trends) 
[Fig. pi a)] . The irreversible nature of the system at large 
wall strain amplitude is illustrated in Fig. [l|c) by peri- 
odically sampled particle tracks over 35 cycles (see also 
Video2 [H]). We note that the largest fluctuations oc- 
cur slightly away from the channel walls, as previously 
observed in steady flows [T7] . 

Irreversible behavior does not occur only where the lo- 
cal strain amplitude surpasses a critical value 7 > 7 C 



[H [5] . As the wall strain amplitude is increased, regions 
in the center of the channel (y/a = 0), where the lo- 
cal strain is negligible, show an increase in diffusivity of 
two orders of magnitude [Fig. [3^a)]. Although the rate 
of increase is greater near the wall, the local diffusivity 
shows a nearly exponential increase with the wall strain 
amplitude including the centerline, where the local strain 
amplitude is negligible [Fig. [3](b)] . 



B. Transition to Irreversibility 

The transition to irreversibility in uniformly sheared 
particle suspensions has been shown to exhibit the char- 
acteristics of a continuous phase transition. The char- 
acteristic relaxation time r for periodically sampled par- 
ticle fluctuations to reach a steady state was shown to 
diverge at the onset of irreversibility in the case of uni- 
formly strained suspensions 2 a . Following the analy- 
sis of Corte et al., the relaxation time is obtained by 
fitting a time series of the local particle MSDs con- 
sisting of 500 cycles from periodically sampled data to 
/ (*) = (/° - f°°) e~ t,T /-t a + f°°, where a = 0.6 and f° 
and /°° are the initial and steady values of the particle 
MSD, respectively [2]. Figure |2Fa) (inset) shows sam- 
ple time series for both low and high strain with fits for 
y/a = 0.44 (see also Video2 [IT]). 

For small wall strains, particles reorganize quickly 
within a few cycles (r < 10 cycles) and settle into a 
reversible steady state. As the wall strain is increased, r 
increases abruptly across the entire channel at 7q m 3.4 
indicating a simultaneous onset of irreversible behavior 
[Fig. Qa)]. Thus, the local critical strain is elevated 
at the wall and depressed near the channel center (com- 
pared to 7 C = 1.6 for 4> = 0.4 in uniform shear) pQ. We 
also note that r does not diverge, but rather maintains 
an elevated value (« 130 cycles) as 70 increases further. 

The transition to irreversibility is sensitive to the par- 
ticle concentration [I], and for the conditions tested here, 
the volume fraction can vary from 30% at the walls to 
50% in the center of the channel as a result of shear- 
induced migration [Fig. W[c)]. The delayed onset of ir- 
reversibility can be explained partially by this effect. In 
Fig. 4lb), we examine the variation in local critical strain 
amplitude as a function of the local particle concentration 
(for the single detected critical wall strain, 7§) and com- 
pare it to the critical strain amplitude for the previously 
studied Couette flow. The much steeper variation here 
compared with the uniformly strained case shows that 
concentration variations alone are insufficient to explain 
the simultaneous criticality across the channel. 



C. Long-range Correlated Fluctuations 

An important additional factor is the role of irre- 
versible activity near the walls, which appears to influ- 
ence the core of the channel some distance away. As 
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FIG. 4: (color online) (a) Particle fluctuation relaxation time 
increases abruptly with wall strain amplitude across the en- 
tire channel. Inset: Sample mean square displacement time 
series for 70 = 2.8 (□) and 70 = 4.0 ( ) at y/a = 0.44 
with fits, (b) The local critical strain, 7 C , observed at the 
global onset 70 ~ 3.4, as a function of the local concentra- 
tion, <j). The much steeper variation here compared with the 
uniformly strained case demonstrates that concentration vari- 
ations alone are insufficient to account for the simultaneous 
criticality across the channel. 



evidence for this effect, we note that coordinated parti- 
cle motion is observed in the center of the channel once 
7o has been exceeded [Fig. RJc)] (see also Video2 [II]). 
Other works have drawn attention to the importance of 
correlated particle motion or "transient clusters" to the 
rheology and flow of suspensions [35] . We quantified this 
collective motion by a displacement correlation function 
g(r') — (Arj • Ar-j)/(A?7 • Ar-j), where r' = \fj — rj| is 
the separation between any two particles. The corre- 
sponding correlation length is denoted by £; the largest 
measurable correlations are limited by the finite system 
size. Below the critical strain, particle motion in the pe- 
riodically sampled data is minimal and spatially uncor- 
rected across the channel [Fig. [5] (inset)]. However, once 
the critical strain is surpassed, the correlation length in- 
creases in the center, but remains small near the walls 
[Fig. [5]. The larger correlation length of the particle 
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FIG. 5: (color online) Particle motions in the center of the 
channel become correlated over a longer range above the crit- 
ical wall strain amplitude, which indicates that the fluctu- 
ations do not arise from local interactions. Inset: Sample 
correlation length profiles for 70 = 2.3 (□) and 70 = 4.0 ( ). 



motion in the center region implies that these fluctua- 
tions originate from irreversible disturbances some dis- 
tance away rather than from local particle interactions. 



IV. CONCLUSION 

We have studied the onset of irreversible motion in 
particle suspensions in an oscillatory channel flow, where 
the shear strain is non-uniform. The transition is shown 
to be quite different from the previously studied Couette 
flow case in that it does not exhibit a divergent relax- 
ation time. Strikingly, the transition is delayed to higher 
wall strain amplitude, and occurs simultaneously across 
the channel, even where the local strain is negligible. We 
find that self-organization and shear-induced migration 
alone cannot explain the simultaneous onset. The tran- 
sition to irreversible dynamics is not governed purely by 
the local strain, but depends also on disturbances gener- 
ated some distance away. This work is an important step 
in understanding the coupled behavior of irreversible par- 
ticle motion (activity), local particle concentration, and 
imposed strain 4 . However, many challenges remain 
including consideration of nonlocal effects in models of 
the transition, the possible role of finite system size, and 
other materials such as polydisperse or soft suspensions. 
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